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Q Abstract: This is a truncated version of the paper Degenerate diffusion with a drift potential: 

a viscosity solutions approach, co-authored with I. C. Kim. The purpose of this version is to 
withdraw the claim of quantitative rate of convergence of the free boundary on the part of 
H. K. Lei. The difference from the previous version lies in Section 3 where 1) the quantita- 
tive version of the convergence of the free boundary statement has been removed and 2) the 
0^ more basic version of some convergence of the free boundary given uniform convergence of the 

function has been rewritten. 

^— > 

We introduce a notion of viscosity solution for a nonlinear degenerate diffusion equation with 
a drift potential. We show that our notion of solution coincides with that of the weak solu- 
tion defined via integration by parts. As an application of the viscosity solutions theory, we 
show that in the case of a strictly convex potential, the free boundary uniformly converges to 
equilibrium as t grows. 

o 

1 Introduction 

Consider a C 2 -potential &(x) : R n — > M, and consider a nonnegative, continuous function 
po(x) : M. n — > M. which has compact support Qq. In this paper we study the porous medium 
. ^ equation with a drift 

X Pt = A(p m ) + V-(pV^, (1.1) 

H 

for m > 1, with initial data po(x). It will be convenient to change from the density variable 
to the pressure variable 

u=—o m ~\ u Q = ^-p™-\ (1.2) 
m — 1 m — 1 

so that the equation becomes 

(PME-D) u t = (m- l)uAu + \Vu\ 2 + Vu • V$ + (m - l)uA$ 
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(for more on the density to pressure transform see e.g., the discussions in [3j). We consider 
continuous and nonnegative solutions in the space-time domain Q = W 1 x (0, T) for some 
T > 0, with prescribed initial conditions u(x, 0) = uq{x) 6 C(R n ). 

When $ = 0, (PME-D) is the widely-studied Porous Medium Equation (PME): We 
refer to the book p3] for references. Moreover, when V = |x| 2 , (PME-D) is obtained as a 
re-scaled version of the (PME) via the transform 

6>(?7,t) := t~ a u(x,t),7] = xt' l3 ,T = hit 

(here u solves (PME)). This suggests that the local behavior of (PME-D) is similar to that 
of (PME), with perturbations due to the inhomogeneity of <3?. We will illustrate this fact in 
the construction of various barriers in Section 2. 

The weak solution theory for (PME-D) in the case of bounded domains has been de- 
veloped in [2] and |6j. Also, in [7], existence and uniqueness of solutions are established 
for the full space case under reasonable assumptions (either the initial data is compactly 
supported or the potential has less than quadratic growth at infinity). 

Further, uniform convergence to equilibrium for (PME-D) has also been shown in [2] (see 
Theorem 3.1). In [8], the connection between the (PME) and the nonlinear Fokker-Planck 
equation is established, which facilitates the use of the entropy method to derive an explicit 
L 1 rate of convergence. In [9], an extensive study is made of a general form of the nonlinear 
Fokker-Planck equation, i.e., pt = V • (Vip(p) + pV^) with suitable assumptions on tp and 
exponential L 1 rate of convergence is obtained. (PME-D) falls under the framework of [9], 
and in fact it is the case that almost all of our results would also go through for a general 
equation of this form, but for ease of exposition we will restrict attention to (PME-D). 

We introduce a notion of viscosity solution for the free boundary problem associated 
with this equation, which we will show to be equivalent to the usual notion of weak solutions 
- see [5] for the general theory of viscosity solutions. Note that, formally, the free boundary 

r(«) := d{u > 0} 

moves with the outward normal velocity 

Ut Vu ._ . _, Vu 

V = — — = {Vu + V$) • t=— = \Vu\ + V$ • 



\Vu\ |Vn| |Vit| 

where the first equality is due to the fact that u = on T(u). In this regard we closely 
follow the framework and arguments set out in [3] (see also [T0j and |3j), where the viscosity 
concept is introduced and studied for the Porous Medium Equation. We point out especially 
that [3j extends the result of [3] to the case where the diffusion term is multiplied by more 
general nonlinearities; our focus, however, is on the added drift term, which introduces 
spatial inhomogeneities. The key utility of the viscosity concept here is that we will be able 
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to describe the pointwise behavior of the free boundary evolution by maximum principle 
arguments with local barriers. As an application, we are able to extend the results of [2] 
and [9] to a stronger notion of free boundary convergence. We summarize our main result 
in the following: 

Main Theorem. There exists a viscosity solution u of (PME-D) in the sense of Definition 



2.1 and 



(a) u is unique and coincides with the weak solution studied in and 12/ . 



(b) Suppose $ is strictly convex, then we have convergence of the free boundary: T(u(-,t)) 
uniformly converges to r(itoo), where Uoo is the unique equilibrium solution to which 



u tends as t — > oo (see Theorem 3.1). 



We will separately (re) state and prove items (a) and (b) in the above as Theorem 2.24 
and Theorem 



Remark. 



o The free boundary convergence may not hold if |V$| vanishes at some points, even 
though the uniform convergence of the solution still holds. This is the reason for 
assuming <E> to be strictly convex. 

o In the case of <S>(x) = \x\ 2 (that is for the renormalized (PME)) Lee and Vazquez |12j 
showed that the interface becomes convex in finite time. It is unknown whether such 
results hold for general convex potentials: we shall investigate this in an upcoming 
work. 



2 Viscosity Solution 

In this section we introduce the appropriate notion of viscosity solution for (PME-D) and 
show that it is equivalent to the usual notion of weak solution. Our definition descends from 
those in [3] and [10]. For more details we also refer the reader to the definitions, discussions 
and results in [3J. 

2.1 Definition and Basic Properties 

Let Q := VL n x (0, oo). For a nonnegative function u(x,t) in Q, we define the positive phase 

n(u) = {u > 0}, tt t (u) ■= {x : u(x,t) > 0} 
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and the free boundary 

r(u) = an(u), r t ( u ) :=dn t ( u ). 

As in [3], to describe the free boundary behavior using comparison arguments we need an 
appropriate class of test functions to handle the degeneracy of (PME-D). 

Let S be a smooth, cylinder-like domain in W 1 x [0, oo), i.e., 

E = (J £(i) x {t}, where £(i) is a smooth domain in M n . (2.1) 

t\<t<t 2 

Definition 2.1. A nonnegative function u G C 2,1 ({ii > 0} n £) is a classical free boundary 
subsolution in £ if 

o n satisfies (PME-D) with < replacing = in the classical sense in {u > 0} n £; 
o |-Du| > on r(n) n E with outward normal velocity 

Vu 

V < |Vu| + V$ • r on Tfu), 
| Vn| 

or, equivalently, 

u t < \Vu\ 2 + V$ • Vn onr(«). 

We define a classical free boundary super solution by replacing < with >. 

Finally, u is a classical free boundary solution if it is both a sub- and supersolution. 

Before proceeding further it is convenient to introduce some auxiliary definitions. 

Definition 2.2. Let (p be a continuous, nonnegative function. Now if ip is another such 
function, then we say that ip touches ip from above at (xo,io) in E if y5 — ^ has a local 
minimum zero at (xo,to) in E n {t < to}- We have a similar definition for tp touching ip 
from below. 

Definition 2.3 (Strictly separated data). For two nonnegative functions u,v : W 1 — > R, 
we write uq ~< vq if the following holds: supp(uo) is compact and supp(uo) C Int(supp(uo)) 
and inside supp(uo), uq(x) < vq(x). 

We note that e.g., due to the maximum principle, a classical free boundary subsolution 
that lies below a classical free-boundary supersolution at time t\ > cannot cross the 
supersolution from below at a later time ti > t\. This observation leads to a notion of 
viscosity solution which takes into account the free boundary: 

Definition 2.4. Let u be a continuous, nonnegative function in Q. 
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o u is a viscosity subsolution of (PME-D) if, for any given smooth domain E given in 
(2.1), for every ip £ C 2,1 (E) that touches u from above at the point (xo,to)> we have 

(p t <(m - l)(pAip + |V(/?| 2 + V<£ • V$ + (m - l)<pA$. (2.2) 

o n is a viscosity supersolution of (PME-D) if, for any given smooth domain S as given 
in (2.1), 

(i) for every <p £ C 2,1 (E) that touches u from below at the point (xo,io) S O(u) D E, 
we have 

<p t >(m- l)ipAp + |V(^| 2 + V99 • V$ + (m - l)<pA$. (2.3) 

(ii) for every classical free-boundary subsolution ip in E, the following is true: If 
i/j -< u on the parabolic boundary of E, then ip < u in E. That is, every classical 
free-boundary subsolution that lies below u at a time ii > cannot cross u at a later 
time ti > t\. 



o u is a viscosity solution of (PME-D) with initial data uo if it is both a super- and 
subsolution and u uniformly converges to u$ as t — > 0. 

Remark 2.5. /n general one can define viscosity sub- and supers olutions respectively as 
upper- and lower semicontinuous functions. Such a definition turns out to be useful when 
one cannot verify continuity of solutions obtained via various limits. This problem does not 
arise in our investigation here thanks to J|^, and therefore our definition assumes continuity 
of solutions. 



It is fairly straightforward to verify that a classical free boundary sub- (super)solution 
is also a viscosity sub- (super)solution. 

Lemma 2.6. If w is a classical free boundary sub- (super) solution to (PME-D), then w 
is also a viscosity sub- (super) solution. 



Proof. We will be brief: The subsolution case presents no difficulty since if contact with 
some ip £ C 2,1 (E) occurs in Q(w) then we use the fact that w is classical there, whereas 
no contact can occur on the free boundary unless |Vto| = 0, in which case the differential 
inequality is satisfied since then <p = |Vy| = and (ft < 0. 



If w is a classical free boundary supersolution, then (i) in Definition 2.4 follows as before. 
To see (ii), let us note that if (p is a classical free boundary subsolution which crosses w, then 
since the free boundary is C 2 , Hopf's Lemma implies that at the touching point | V(/?| < \Vw\ 
(see e.g., p2]). On the other hand, since (p started below w, at the touching point we must 
have v n (ip) > v n (w), which leads to a contradiction since it is also the case that we have 



|Vu>| \Vip\ ' 



□ 
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Next we have the following stability result. 

Lemma 2.7. Let u £ be a smooth solution of (PME-D) with initial data uq + e and let u be 
its uniform limit. Then u is a viscosity solution of (PME-D) with initial data uq. 



Proof. Let E be as given in (2.1) and let tp £ C 2,i (S). 



1. Let us first show that u is a subsolution. First suppose that (p touches u from 
above at the point (xo,to). We may assume that u — ip has a strict maximum at (xo,io) i n 
Q(u) n S n {t < to} by replacing <p by 

<p(x, t) := ip(x, t) + a((x - x ) 4 - (t — t ) 2 ), a > 

if necessary. By uniform convergence there exists a sequence (x £ ,t £ ) converging to (xo>*o) 
such that u £ — p has a local maximum at (x e ,t e ). Now if we we let 

ip(x, t) := p(x, t) - ip(x £ ,t £ ) + u e (x £ , t £ ) 

Then u e — p has a local maximum at (x £ ,t £ ) with (u £ — (p){x e , t £ ) = 0. We can now conclude 
by taking the limit of the viscosity subsolution property of u £ . 

2. Next we show that u is a supersolution. Let (p be a classical free-boundary subsolution 
such that <p(x, ii) -< u(x, t\). Since the u^s are strictly ordered, u < u £ and hence <p(x, t\) ~< 
u(x,ti) < u £ (x,ti). Now suppose p touches u £ at some point (x2,t2), then p(x2,t2) > 
since u £ is positive, so by continuity, there is a parabolic neighborhood of (^2,^2) in which 
both functions are classical and positive. By the Strong Maximum Principle, the touching 
cannot have occurred at (#2, £2)5 a contradiction. We conclude that p> < u £ so that in the 
limit p> < u £ . 

□ 

An immediate consequence of the above lemma is that weak solutions are viscosity 



solutions (see Corollary 2.12). We shall introduce the precise notion of weak solutions in 



the next subsection, and summarize some results from [2]. 
2.2 Weak Solutions 

To be consistent with the setup in both [2] and [9] , let us return to the density variable and 



consider the solution of (1.1) in a bounded domain Q with Neumann boundary condition: 



(AT) 



Pt = Ap m + V • OV$) in ftxR+, 
dp m /dv + p (d$/dv) = on dQ x R+, 
k p(x,0) = po(x) in 17. 
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We will see shortly that we need not worry about the fact that we are on a bounded domain, 
but for now we will let Q = £1 x R + and Qt = Q x (0, t]. As in [2], we make the following 
definition: 

Definition 2.8. We say p : [0, oo) — > L 1 (0) is a weak solution of (N) if 

(i) p G C([0,t];L 1 (fi))nL oo (Q i ) for all t € (0,oo); 

(ii) for all test functions ip G C 2,1 {Q) such that ip > in Q and dip/dv = on dQ x R + , 
we have 

/ p(t)<p(t) = [ p(0)<p(0) + [ [ (p<pt + p m A<p - • 
Jn Jn J jQ t 

We also define a weak subsolution (respectively supersolution) by (i) and (ii) with equality 
replaced by < (respectively >). 

From pj we have existence, regularity, uniqueness and comparison principle for weak 
solutions: 

Theorem 2.9 (From |2j). Under the assumption that <3? is C 2 in Cl: 

(a) the problem (N) has a unique solution; 

(b) the solution is uniformly bounded in Q and is continuous in any set £1 x [0,T]; 

(c) suppose p{t) is a subsolution and p(t) is a supersolution, then if p Q < p in f2, then 
p(t) < p{t) innfort>0. 



The existence of solutions is obtained as the uniform limit of solutions to uniformly 
parabolic problems (equicontinuity is obtained from [6]). For our purposes, a very simple 
approximation basically suffices and we summarize the relevant result in the following: 

Lemma 2.10 (From |2j). Let u e be a solution of (N) with initial data Uq = uq + e, then 
u £ is equicontinuous and there exists a subsequence which uniformly converge to u which is 
the unique weak solution to (N) with initial data uq. 



While a priori our viscosity solution is defined in all of W 1 , since (formally at least) 
solutions of (PME-D) should have finite speed of propagation, the boundary conditions 
should be inconsequential if we take sufficiently large. (Later we will also establish finite 
propagation for viscosity solutions - see Corollary 2.16 ) Control on the speed of expansion 
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of the support can be done via comparison with any (weak) supersolution. In particular, 
when <3? is convex, we can use the stationary profiles of the form ^(x) = (C — $)+ with 



sufficiently large C as a supersolution (see Theorem 3.1) 



Remark 2.11. Alternatively (and perhaps this is a cleaner line of reasoning), we can 
directly use the result of /7/ on existence and uniqueness of solutions in all of W n , which 
implies in particular that the results of f^jj also apply in that setting. 



Combining Lemma 2.7 with Lemma 2.10 and the uniqueness statement in Theorem 2.9 
we obtain: 

Corollary 2.12. Any weak solution is also a viscosity solution. 

We will eventually establish uniqueness of viscosity solutions via maximum-principle 
type arguments, which culminates in the identification of the two notions of solution. 



2.3 Construction of test functions 

In this subsection we collect some test functions, i.e., (classical free boundary) sub- (super) 
solutions, to (PME-D) which will be useful for comparison purposes. In the first couple 



of lemmas (Lemmas 2.15 and 2.18) the idea is to control the <3? dependence via Taylor 
expansion in a small neighborhood of a point, so that we can appropriately perturb the test 
functions for (PME) constructed in [3] and [1] for our purposes. 

Let a > be small. We will first explain how locally the (PME-D) can be viewed as a 
perturbation of (PME). The starting point is to observe that if u is a solution of (PME-D) 
in the parabolic cylinder 

Q a (xo,t ) = B(x ,a) x [t , to + a}. 
By Taylor's theorem, u satisfies an equation 

ut = (m- l)uAu + \Vu\ 2 + b-Vu + (m- l)itA$ + 0(a)[D 2 $]Vu + 0(a 2 ), 

where 

b = V$(x ,t ) 

and [D 2 $] denotes the relevant term in Taylor's expansion for V3>. 
On the other hand, if we use hyperbolic scaling and define 

v(x, t) = a~ 1 u(a(x — xo),a(t — to)), 

then v satisfies 

v t = (m - l)vAv + \Vv\ 2 + Vv ■ V$ + a{m - l)uA<J> 
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in the parabolic cylinder 

Qi (0,0) = 5(0,1) x [0,1]. 
Combining this with the Taylor expansion, we see that v satisfies 

v t = (m - l)vAv + \Vv\ 2 + b-Vu + a(m - l)vA$ + 0(a)[D 2 $]Vv + 0(a 2 ) 
:= (m - l)vAv + \X7v\ 2 + b-Vu + R(x ,v). 

From Taylor's theorem, we see that 

\R{x ,v)\ <a||$||o2(H + |Vi;|). (2.4) 

Finally, we define 

w(x, t) = v(x — bt, t). 

Then uu satisfies 

wt = (m — l)wAw + \Vw\ 2 + R(x , w), (2.5) 
which is indeed the (PME) with an 0(a) perturbation in the unit scale parabolic cylinder 

Qi(bt,0) = B(bt,l) x [0,1]. 

To construct test functions for (PME-D) from perturbations of solutions of (PME) we 
shall have to reverse the order of operations: We will construct super- (sub)solutions of 



(2.5) by perturbing corresponding super- (sub)solutions of the (PME) at the unit scale and 
then we will use hyperbolic scaling to arrive at a rescaled (and translated) local super- 
(sub)solution for the (PME-D). Let us note for future reference that 

u(x, t) = aw (^a^ 1 (x + bt), a -1 i^ . 
The next proposition gives one way to perform the necessary perturbation on solutions 



of (PME) to arrive at (2.5). 



Proposition 2.13. Let u(x,t) be a viscosity subsolution of (PME) in Bi+ a (0) x [—1,1]. 
Then for < a < 1 , 

ui(x, t) := e~ at sup u(y,t) 

y<=B a - at (x) 

is a subsolution of 

(PME-sub) ( u i)t = (m — l)u\Aui + |Viii| 2 — a|Vtti| — au\ 

in -Bi(O) x [—1,1]. Similarly, if u(x,t) is a viscosity supersolution of (PME) in -Bi(O) x 
[-1, 1], then for0<a<l, 



U2{x,t) := e at inf u(y,t) 

y£B a - at (x) 
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is a supper solution of 

(PME-super) i u 2)t = (m — 1)u2Au2 + |V«2| 2 + a|Vu2| + Q.U2 

in Bi(0) x [-1,1]. 

Proof. We only show the supersolution part. Let v>2 be as given above. Suppose then that 
(p is classical and touches U2 from below at some point (xo,to). We first note that there 
exists x\ G i? Q _ a t(xo) such that U2(xo,to) = e atl u(x\, to). 

Next for any unit vector b, let us consider 

<p(x, t) = e~ at ip(x — (x\ — xo) — ab(t — to),t). 

Then we note that 1) <p(x±,to) = e~ at °ip(xo,to) and so (u — <p)(xx,to) = and 2) by the 
definition of 112 as an infimum and by continuity of ip, in a small parabolic neighborhood of 
(xo, to), it is the case that u — (p > U2 — (p > 0; we therefore conclude that (p touches u from 
below at {x\,to) and so we have, 

[<pAp + \V<p\ 2 }(x ,to) = e at »[^A<^+|V^| 2 ](a;i,to) 

< e ato Mxi,to) 

= [<p t ~ oup - ab -V(p](x ,to). 

Now the desired inequality is achieved by setting b = |^r(^0)*o)- 

Indeed the above calculation shows that if u is a viscosity supersolution of (PME), then 
U2 should be a viscosity supersolution of (PME-super): If a classical free boundary subsolu- 
tion tp of (PME-super) crosses 112 from below, then the corresponding (p is a subsolution of 
(PME) and crosses u, yielding a contradiction (there is no distinction between the interior 
and boundary cases). □ 

We will use the spherically symmetric supersolutions for (PME) from [3]: 
Lemma 2.14. Consider the function 

H(x,t;A,oj) = A(\x\ + ojt - B)+ 

with R/2 < B < R. Then u is a classical free boundary supersolution of (PME) in the 
domain {\x\ < R} x [^{B - R),0] if 

to , w ,R-B 

_>l + 2(m-l)(n-l)— — . 
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Proposition 2.13 and Lemma 2.14 yield the following: 



Corollary 2.15. Let us fix xq £ M n and let H be given as in Lemma 2.1^, Then the inf 
convolution of H , given as 



H{x,t;a) = e at inf H(y,t) 

y£B a - a t(x) 

is a classical (free boundary) supersolution of (PME-super). Consequently, there exists 
C = Cq which only depends on the C 2 -norm of & in B\(xq) such that 

H(x, t) := aH_(a' 1 (x — xq + b(t — to)), a~ l {t — to); Ca) 

is a classical (free boundary) supersolution of (PME-D) in Q a := B a {xo) x [to — a, to]. 



Proof. By Lemma 2.6 and Proposition 2.13, H_ is a viscosity supersolution of (PME-super) 
and hence adjusting a by multiplying by a suitable constant to take into account the C 2 - 



norm of we obtain a viscosity supersolution of (2.5 ) from which we obtain a supersolution 
of (PME-D) in Q a . 

To finish it is sufficient to show that H has the required regularity to be a classical free 
boundary supersolution. For this, note that for simplicity we have only taken the supremum 
over space and the reader can readily check that in this case the infimum for H_ (x, t; a) is 
achieved at the point y which minimizes \y\ subject to the constraint that \y— x\ = a — at, 
and thus an explicit expression is possible for H_ . □ 



By comparison with these supersolutions, we immediately obtain 

Corollary 2.16. Any viscosity solution has finite speed of propagation and is bounded in 
a big ball in any local time interval. Further, if $ is convex, then via comparison with the 
stationary solutions of the form given in Theorem\3.1\ the above holds globally in time. 



For subsolutions we will make use of the Barenblatt profiles (see e.g., [3] and |14j). 
Lemma 2.17. [Barenblatt] Let B{x,t;r,C) be the family of functions 

B( ,,«,T. C ) = < c <' + f -*l'''>+ 

[t + T) 

with constants X,K,C,t > such that 

\ = ((m - \)d + 2)" 1 , 2K = A. 
Then B(x,t;r,C) is a classical (free boundary) solution of (PME). 
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Using Proposition 2.13 (see also the proof of Corollary 2.15) once again, we obtain the 
following: 

Corollary 2.18. Let us fix some (xo,to) and let B be a Barenblatt function. Then there 
exists a constant C which only depends on the C 2 -norm of & in B\(xq) such that 



B{x,t) 



ae 



-Ca(t-to) 



sup B(a (y 



x 



a^(t-to)),a -1 (*-*o)) 



is a classical (free boundary) subsolution of (PME-D) in Q a := B a (xo) x [to — en, to]. 
Remark 2.19. The reason for taking the hyperbolic scaling is because we will have occasion 



to require rather fine control on the boundary velocity (see Lemma 2.23) and this is the 
scaling which preserves the velocity - in contrast to the parabolic scaling, which dramatically 
reduces the effect of the drift $ in the bulk (the positivity set), but unfortunately at the cost 
of severely disrupting the boundary velocity. 



To establish the Comparison Principle, we will need the following weak analogue of (ii) 
in the definition of viscosity supersolutions for subsolutions, the proof of which utilizes an 
approximation lemma from [3]. 

Lemma 2.20. Let u be a viscosity subsolution of {PME-D), and let ip be a classical free 



boundary supersolution from Lemma 2.15 which lies above u at some time to- Then ip 
cannot cross u from above at a later time t > to- 



Proof. Let ip and u be as described in the statement, and suppose that tp touches u from 



above at some point (xo,to)- From Lemma 2.15 we have that (p is given as the inf convo- 
lution of some spherical traveling waves from Lemma 2.14[ which we denote by ip. Further, 
let us suppose the infimum is achieve at (xi,ti) so that 1) (p(xo,to) = ip{x\,ti) and 2) by 
the definition of (p as an inf convolution, the translated function 

tp(x, t) = ip(x + (x 1 - X ), t + (ti - t )) 

also touches u from above at the point (xo,to)- From Lemma 4.4 in [3], we know that 
ip can be given as the monotone limit of classical positive supersolutions, and hence the 
same is true of ip: I.e., there exsits ip e \ ip, with %[) e > classical. But since u cannot 
touch ip £ by the Strong Maximum Principle, we obtain in the limit that u < i/j, which is a 
contradiction. □ 



2.4 Comparison Principle and Identification with Weak Solution 

Here the outline of the proof closely follows that of the corresponding result for (PME) 
(Theorem 2.1 in jlj): We will give an abridged version of the proof, pointing out main steps 
and modifications for our problem. Here we let Q = M n x M + . 
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Theorem 2.21. [Comparison Principle] If u is a viscosity subsolution and v is a viscosity 



supersolution in the sense of Definition 2.4 with strictly separated initial data, uq -< vq, then 
u(x,t) < v(x,t) for every (x,t) G Q. 

Proof. 1) [Sup and Inf Functions] For given 5 > and r > small with r>iwe introduce 
the regularized functions 

W(x,t) = _ inf v(y,r) 

B r _ 6t (x,t) 

and 

Z(x,t) = sup u(y,r) 
B r (x,t) 

First note that W and Z preserve properties of v and u: 

o W is a supersolution and Z is a subsolution; 
o Z(-,r) -< W(-,r) for r sufficiently small. 



For a proof of the first item see [3], Lemma 7.1 or the proof of Lemma 2.15 The proof of 
the second item relies on the fact that 

o The support of viscosity subsolutions and supers olutions evolve in a continuous way. 
Here continuity is understood as continuity in the Hausdorff distance (in time) of the 
positivity set. 

The proof of this can be done by comparison with the supersolutions (respectively subsolu- 
tions) constructed in Corollary 2 . 1 5| (respectively Corollary 2.18). We omit the details since 



with replacement of barriers it is no different from the proof of Proposition 6.2 in [3]. 

Thus if we can prove that W stays above Z for all choices of r and 5 (sufficiently small), 
then we may take 5 — > and then r — > to recover the conclusion for u and v. First 
let us note that due to the Strong Maximum Principle, W cannot touch Z from above, 
and therefore we are reduced to the analysis of a first contact point of W and Z at some 
P = (A,t ). 

The key usefulness of Z and W lies in the fact that they enjoy interior/exterior ball 
properties: 

o The positivity set of Z has the interior ball property with radius r at every point of its 
boundary and at the points of the boundary of the support of u where these balls are 
centered we have an exterior ball; 
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Figure 1: The geometry at the contact point 



o The positivity set ofW has the exterior ball property with radius less than r — 5t (since 
in this case we really have an exterior ellipsoid in space-time) and at the points of the 
boundary of the support of v where these balls are centered we have an interior ball. 

For detailed proofs of these statements we again refer the reader to [3J. 

2) [The Contact Point] The first contact point Po = (xQ,to) is located at the free 
boundary of both functions. Therefore by the definitions of Z and W, there is a point 
Pi = H) on t ne free boundary of u located at distance r from Po an d there is another 
point P2 = (x2,t2) on the free boundary of v at distance ro = r — 5to from Po. Let us 
also denote by Hz (respectively H\y) the tangent hyperplane to the free boundary of Z 
(respectively W) at Po. (see Figure [I]) 

Lemma 2.22. Neither Hz nor H\y is horizontal. In particular, one can denote the space- 
time normal vector to Hz, in the direction of P\ — Pq, as (e n ,m) 6 M n x M where \e n \ = 1 
and —00 < m < 00. 

Proof. It is enough to show that t\ > t$ — r (i.e., T(Z) cannot propagate with infinite speed) 
and t2 < to + r (i.e., r(W) cannot propagate with negative infinite speed). The desired 
conclusion then follows by the ordering of Z and W. 
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i = t 2 



P(V) 



lb • 



v > o 



Figure 2: Infinite, negative speed 



We first show that £2 < to + r. Otherwise Hyy is horizontal and after translation we 
have Po = (0; ~ r ) an d P2 = (0,0). Moreover £l(v) has an interior ball at P2 with horizontal 
tangency with radius < r' < r. Now in any parabolic cylinder 

C v = {(x, t) : \x\ <r],-r] 2 <t< 0} 

with bottom edge contained in the interior ball (which can be achieved by taking ij < r'), 
we have by continuity that v > M > on that edge, (see Figure [2]) 



0, so we will have to compare 

2 and 



On the lateral boundary of C v it may be the case that v 
with a subsolution with support strictly contained in —rj < \x\ < ij at time t = —rj 
still contains in its support at time t = 0, which rigorously implies that v cannot contract 
sufficiently fast for (0,0) to be a free boundary point. The necessary subsolution can be 



constructed as the one in Lemma 2.18, adjusted for the parabolic scaling. 



The case ti > to — r follows mutatis mutantis from the arguments in [1] and [3] (see [1], 



Lemma 4.2 or [3], Lemma 8.2) with the barrier in Lemma 2.15 replacing the barriers used 
therein. □ 



3) [Non-tangential Estimate] The next lemma states that the normal velocity V of T(Z) 
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at (xQ,to) satisfies, in the viscosity sense, 

V ixo ,t ) < (\VZ\ + V$(x ) • ^\ (x ,t ). 

Lemma 2.23. Let x n '. — x • e n , and consider a non-tangential cone 

K := {x : x n > k\x\ with k > 0}. 

Then we have 

Z(x + x,t ) 
nm > m — \/<p{xq) ■ e n . 

X&K, x^O x n 

Proof. The argument is parallel to the proof of Lemma 4.3 in [3]; the only difference for 
us is taking into account the change of reference frame introduced by the drift given by 
This is ensured by the local nature of the construction of our barrier in Lemma 2.15[ which 
replaces the corresponding barriers used in [3]. 

□ 



4) [Conclusion] Due to Lemma 2.23, we may place a small subsolution (p from Lemma 



2.18 below Z at Pq with speed close to m (again see Remark 2.19, which assures us that 



our subsolutions are constructed so that this is possible) such that it crosses anything with 
speed m! < m. Since (p is also below W and hence v (after a small translation), v must 
expand by at least m' , but then T(W) has speed m' + 5 > m at Pq, yielding a contradiction 
to the fact that Z touched W from below at Pq- D 



We can now establish uniqueness of viscosity solutions: 

Theorem 2.24. The problem (PME-D) admits a unique solution in the class of viscosity 
solutions as defined in Definition \2.J\ for continuous and nonnegative initial data. This 
solution coincides with the continuous weak solution. 



Proof. The existence of a continuous weak solution can be provided as the uniform limit of 



classical solutions with initial data uq :£ = uq + e, and by Lemma 2.7 such a limit, which 
we will denote by U, is also a continuous viscosity solution. Further, by comparison with 
UQ tE and taking a limit, it is clear that such a limit U is also a maximal viscosity solution. 

Uniqueness would follow if we can show that any other viscosity solution u also cannot 
be smaller than U . For this purpose, consider u n (x, t) with initial data u n (x, 0) := (uo — ^) + . 
Now consider positive u^ 1 such that |u^ n — u n \ < - in M. n x [0, T]. It follows from Lemma 



2.10 that u^ 1 uniformly converges to U2(x,t), which is then a continuous weak solution of 



(PME-D). Therefore, by uniqueness of weak solutions, U2 is equal to U. On the other hand 



by Theorem 2.21 u n -< u and thus U = U2 < u. Hence we conclude. □ 
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Using Theorem 2.21 and Theorem 2.24, we can in fact prove a stronger comparison 
theorem for viscosity solutions (see [3], Theorem 10.2): 



Theorem 2.25. Let u\ and U2 be respectively a viscosity subsolution and a viscosity super- 
solution of (PME-D) with initial data no i and uq^ such that uo,i < ^o,2- Then u\(x,t) < 
u 2 (x,t). 



We can now strength Lemma 2.20 



Lemma 2.26. Let u be a viscosity subsolution of (PME-D), and let ip be a classical free 
boundary supersolution which lies above u at some time to- Then ip cannot cross u from 
above at a later time t > to- 



Proof. Let us first replace <p by a viscosity solution of (PME-D) with the same initial data 
which we denote v. Since ip is a viscosity supersolution by Lemma 2.6 we have that <p > v 
by Theorem 2.25 Finally, u < v by Theorem 2.21 □ 



3 Convergence to Equilibrium 

We begin by discussing the set of equilibrium solutions to (PME-D) and reviewing some 



known results. Since by Theorem 2.24, the unique viscosity solution coincides with the 
continuous weak solution, we may invoke the results of [2] which are stated in terms of 
weak solutions. 

The set of equilibrium solutions and uniform convergence of solutions to the equilibrium 
are established in [2]. We state below a restricted version of the result for our purposes. 

Theorem 3.1. [Theorem 5.1, \2i] Functions of the form (C — $)+ with C£R are equilib- 
rium solutions for (PME-D), and given uq, there exists a unique > such that u(x,t) 
uniformly converges to (Coo — &)+ as t — )• oo. 

It is fairly immediate that the set {(C— $)+, C G M} is contained in the set of equilibrium 
solutions; the converse containment and the convergence statement are established based on 
a ^-contraction result in [2] . The uniqueness of the equilibrium solution comes from the 



fact that due to its divergence form, the density function p(-,t) given in (1.2) preserves its 



L 1 -norm over time. It follows that the equilibrium solution is determined by the condition 

[u } 1/(m - 1) (x)dx = [ [u OQ ] 1 ^ m ~ l \x)dx. 



We can now state our free boundary convergence theorem. 
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Theorem 3.2. Let u OQ (x) be the unique equilibrium solution given by Theorem 3.1 Then 
Tt(u) uniformly converges to r(uoo) in the Hausdorff distance, as t — > oo. 

Proof. It will be convenient to look at the sublevel sets 

${< r } := {x : <£(x) < r}. 

We note that since <I> is strictly convex, the sublevel sets with r £ M start from a point 
(where the global minimum of $ is achieved) and exhaust M n (since |V$| > except at the 
minimum). Let e > 0. Invoking continuity of <3? if necessary, it is sufficient to show that 
there exists a T > such that 

for all t > T. 

Ft{u) C [ < ^ > {<c 00 -e}] C: This is the statement that the free boundary converges from the 
inside and is immediate from the uniform convergence to u^. Indeed, uniform convergence 
implies that e.g., there exists T > such that u(x, t) > e/2 for all (x, t) € ^ > {<c 00 -e} x ITj 
it necessarily follows then that T t (u) lies outside ^{<c x ~e} since otherwise by continuity 
we can find a sufficiently small ball inside < ^{<c 00 -e} where u(x,t) < e/2, contradicting the 
definition of T. 

r t (u) C ^{<c 00 +e} : I n principle this should follow from a "dual" sort of argument, i.e., 
some statement to the effect that if y £ Tt(u), then there exists some sufficiently small 
5 > such that u(-,t) is "large" in Bg(y). (E.g., an analogue of Corollary 2.2 in [I] for 
(PME-D) would suffice for us.) We will instead directly make a comparison argument to 
show that the potential appropriately "pulls" the free boundary inwards, taking advantage 



of the explicit form of the supersolutions constructed in Corollary 2.15 



• The uniform convergence of u(-,t) to (Coo — &)+ means that there exists some con- 
tinuous decreasing function n(t) — > such that 

u(x,t) <n(t), for x (£ ${< Cco }- 

Let us define 

r* = sup{t:J2 t C*< Coo+ ^( t )}, (3.1) 

where fj is a smoothed version of n so that it is differentiable (obtained e.g., by 
convolving with a small smooth bump function) and K > is some (large) constant. 

Let us first make a heuristic, non-rigorous argument to the effect that T* = oo. Let 

x* G d<5>{< Coo+ r(T*)} nr T .(u); 
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such an x* exists by the minimality of r(T*), and morally it represents the worst case 
scenario: x* is the furthest away from where it should be among all points x £ Tt*(u). We 
note that V<3?(x*) and Vit(x*) point in the outward and inward normal directions of IV* (u), 
respectively, and hence have opposite sign. We conclude then that the normal velocity of 
the boundary at x* is 

V(x m ) = |Vu(x*)| + V$ • * . = |Vu(x*)| - |V$(x*)|. 

If we assume for the moment that |Vu(x*)| « u(x*) < rj(T*), and that |V$(x*)| > A, then 
we see that 

V(x*)< V (T*)-A. 



Now if it can be arranged so that 



T,{T*)-A<T?(t), 



then we can conclude that the most outlying point of Tj-* is moving inwards at a faster 
rate than indicated by the rate defining T* , thus yielding a moral contradiction to the 
maximality of T* . 



to 



A rigorous version of the argument will be established by invoking Lemma 2.26 
compare u with a (classical free boundary) supersolution constructed in Corollary 2.15 
and thus all considerations regarding the size of u, Vu and the boundary velocity will be 
transferred to the relevant supersolution. These supersolutions are rescaled and translated 
versions of spherically symmetric supersolutions of (PME) supported in annuli. We remark 



that since we are not establishing a precise quantitative rate, the precise form given in (3.1 ) 
is only for convenience: rj(t) sets a convenient small scale. 



• We will now take xo € dQi^c^+KfjfT*)} to be arbitrary. Let us first describe the geo- 
metric setting in which we will make our comparison argument. Since d^^^^^j^fj^p*^ 
is convex, there is an exterior ball at xo which is centered at some xi, which we denote 
B a (x\). The comparison argument takes place in the parabolic cylinder 

S=[ J B Q (x 1 )n J B a /(x )] x [T*,T* + t]. 

(See Figure[3j) Here a' is sufficiently small so that u(x,T*) < rj(T*) in B a i(xo) and a 
is correspondingly small, and, finally, r is sufficiently small in a way to be computed 
shortly. 

• By translation, we will assume without loss of generality that the relevant barriers 
are centered at (0,0) (of course, really, to = T* , and we shall recover T* when it is 
relevant) and we will denote by e.g., $ = $(xo), with xo as before and taken to be 
the point at which the relevant Taylor expansion is performed. Now returning to the 
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proof of Corollary 2.15, we see that solving the relevant minimization problem gives 



the explicit expression 

H{x, t; Co) = Ae Cat (\x - (Cq(1 - t))x\ + wt - B)+, 

where we recall that C is a constant depending on the C 2 -norm of and B is a 
number of order unity. Thus, 

H(x, t) = aH(a~ 1 (x + bt), a _1 t; Ca) 

= Ae Cat (\x - Ca(l - t)x + bt\ +cot- Ba)+ 

= Ae Cat ( \(\x\-Ca(l-t))x + bt\+ut- Ba )+, 

where b = V<£. So we need to ensure that 



H(x, t) > Ae Cat { [\x\ + ut] - [Ca(l - t) + \b\t + Ba] 
on the parabolic boundary d"E. 

Next we compute the free boundary velocity of H. We have 

H t = CaH + Ae Cat (u> + ^|(|x| - Ca(l - t))x + bt\ 

= CaH + Ae Cat (oj + (y + bt, aty + b^J , 

where y = (\x\ — Ca{\ — t))xisa shrunken version of x, whereas 

Vi7 = V|y + bt\ = y + bt 
has absolute value equal to unity. Thus 
Ht 



V(H) 



\VH\ 

Ae Cat ( Ca{\y + bt\ + uit - Ba) + co + ( y + bt, aty + b 



(3.2) 



= Ae Cat ( (1 + Cat)u> + Ca\y + bt\ + M+ ? 2 (y, 6) + _ 5Ca 2 ] 

(3.3) 

To ensure that locally the free boundary is shrinking faster than the rate given in the 
definition of T* we need that 

V{H) < fj'(T*). 
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We will now amalgamate our conditions and observe the correct choice of parameters. 
Let us first uniformly bound 

|V<D|>A in [<% Coo }] C - 
We will take t < r = 0(a) and note that e.g., 

\y + bt\ > \\x\ - Ca(l - t) - 5t\ ~ ||x| - (C + 6)a\ +0(a 2 ), 



so the expression (3.3) becomes (with A = 1) 

V(H) = uj- c(y, b) + 0{a 2 ) > u - cA + 0(a 2 ), 

for some constant c > of order unity. Here the negativity of (y, b) comes from the 
fact that V$ is pointing outwards, whereas the vector from x\ to x$ is pointing in 
the opposite direction (again see Figure [3]). If we also take fj'(T*) to be higher order 
in a, then we arrive at the condition 

u - cA < 0. (3.4) 



Similarly, taking r/(T*) also to be higher order in a, the condition (3.2) becomes 

\x\+oja> (C + B + 5)a + 0(a 2 ), (3.5) 
where S is some local upper bound for |V$|: |6| < 5. 



Writing \x\ = Xa, the conditions (3.5) and ( |3.4[ ) finally become (for a sufficiently small 
and again with A = 1) 

A + oj > C + B + 5 and cA > u, 

which can be satisfied by taking A to be a large enough factor (this corresponds to 
inserting a factor A in the relevant hyperbolic (re)scaling). 

Since d& {<Coo+Kri(T*)} is compact, we can take a to be uniform over all relevant 
xq £ ^ < ^{<c 00 +-f^»?(T*)} (recall that a basically needs to be small enough to ensure 
that u(x,T*) < r]{T*) in a ball of size a' ~ a) thus yielding a contradiction to the 
definition of T* . 

We now state the result positively. The argument in fact shows that if T satisfies the 
condition 

sf(T) : n T c * { < Coo +Jc«(r)}. 

then t satisfies ^(t) for all t >T. Now by comparison against an equilibrium solution 
of the form given in Theorem 3.1 with a sufficiently large C, it is the case that we 
can choose K > large enough so that there is some time T satisfying £f (T). We can 
now conclude by choosing T even larger (if necessary) so that Kfj(T) < e. 



□ 
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